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The wave-vector integrated dynamical spin susceptibility X2d{<^) of 
YBa2 Cus Oe+y cuprates is considered. X2D is calculated in the superconduct- 
ing state from a renormalized mean-field theory of the t-t' -J -model, based 
on the slave-boson formulation. Besides the well-known "4 1 meV resonance" 
a second, much broader peak ('hump ') appears in lmx2D ■ It is caused by 
particle-hole excitations across the maximum gap A.^ . In contrast to the 
resonance, which moves to lower energies when the hole filling is reduced 
from optimal doping, the position of this 'hump' at ~ 2A^ stays almost 
unchanged. The results are in reasonable agreement with inelastic neutron- 
scattering experiments. 

PACS numbers: 71.10.Fd, 74.25. Ha, 74.72.Bk, 75.20.Hr 



1. INTRODUCTION 

The most prominent feature in the magnetic excitation spectrum of 
YBaaCusOe+y (YBCO) and RjoSrgCaCuaOs+s (BSCCO) cuprates is the 
so-caUed "41 meV resonance"liJa'Ba'Ela at the antiferromagnetic (AF) wave 
vector q = (vr, vr) . Its energy u;p»<, is « 40 meV in optimally doped samples 
and decreases with underdoping3iJ§0 down to cj^es ~ 24 meV . Recently the 
magnetic response has also been studied by averaging the neutron-scattering 
data over the in-plane 2D Brillouin zoneEHill. The resulting local magnetic 
excitation spectrum Imx2D('^) shows the above-mentioned resonance and 
a second, hump-like feature at an energy ujhump above LOres ■ In contrast to 
iOres , ^hump depends only weakly on the doping level. Within the calculation 
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to be presented in the following the 'hump' is naturally explained by particle- 
hole (ph) excitations across the maximum d-wave gap A*^ . The energy 
^hump ~ SA'^ comes out almost independent of doping. The resonance, on 
the other hand, emerges from a ph-bound state in the magnetic (spin-flip) 
channel and shows a strong doping dependence. 



2. MODEL AND MEAN-FIELD THEORY 

Our starting point is the doped Mott insulator. We study the J- model 
on a simple square lattice of Cu-3d orbitals for each of the two coupled Cu02 
layers (planes) in YBCO or BSCCO : 

In the subspace with no doubly occupied orbitals, the electron operator 
on a Cu-lattice site u is denoted Cy^ with spin index cr = ±1 ; Sj^ is the 
spin-density operator. A Cu-site is specified through u = [i, I] , where i = 
1 . . . Nl indicates the Cu-position within one Cu02-plane and I = 1,2 selects 
the layer, t^,^' denotes the effective intra- and inter-layer Cu-Cu-hopping 
matrix elements, and Ji,u' the antiferromagnetic super exchange . To deal 
with the constraint of no double occupancy, the standard auxiliary-particle 
formulation Ci^o- = blf^a- is used. The fermion fl^ creates a singly occupied 
site (with spin a), the "slave" boson 6j, an empty one out of the (unphysical) 
vacuum bu\0) = fua\0) = 0. The constraint now takes the form = 
btbu + J2cT flafvcj = 1 • In mean-field theory the constraint is relaxed to its 
thermal average {Qu) = 1 • Together with the number x of doped holes per 
Cu-site, it fixes the fermion and boson densities to 

{l-x) = Y.{flU,) , x = {b%). (2) 

(J 

The derivation of mean-field equations is presented in Ref. |l^ . The 
dynamical spin susceptibility is given in units of ((7/1^)^ as 

X(q,fe,w) = x^(q,'^)cos2(^g2) +x'{q,t^)sin^{^qz) , 

where q is the in-plane wave vector, d denotes the distance of Cu02 planes 
within a double-layer sandwich. For the even (-I-) and odd (— ) mode sus- 
ceptibilities a RPA-type expression is obtained, 

X q,w — , . , , (3) 
i + ^/^qkrH 



Local Magnetic Excitation Spectrum of YBa2Cu306+ 



The irreducible part x**^^ consists of a particle-hole (ph) bubble of fermions 
as is known from BCS theory, 



k,p^ s,s'=±l 



,ee' + AA' 

1 + ss 



EE' 



fjs'E') - f{sE) 
uj + sE - s'E' + ^ ' 



with p = {q_,Pz) and Pz = {0,7r} for the modes {+, — } . Boson excitations 
do not enter x oii mean-field level. Fermions obey an effective dispersion 
£ = e(k,p^) , e' = £(k + q,pz +Pz) , 

e(k,pz) = -2t[cos{kx) + cos{ky)] - At' cos{kx) cos{ky) - t-^Ck^e"^" 

and d-wave gap function A = A(k, pz) , A' = A(k + q^,pz + Pz) , 

A{k,pz) = ^[cos{k^) - cos{ky)] + A^^e'P' 



These enter the usual quasi-particle energy E = Ve^ + A^ , E' = V e'^ + A'^ . 
Formally, vertex corrections to the simple bubble Eq. (Q) have to be taken 
into account. However, these have almost no effect in the energy-range below 
2A° and are therefore ignored in the followingEj . 

From Fey n man's variation principle the effective hopping parameters 
are determined adiHl t ps x t + 0.15 J , t' = xt' , t^(k) « x t"'"(k) . For the bare 
nearest and next-nearest neighbor hopping we assume t = 2J , t' = — 0.45t , 
and for the inter-plane hoppingtJll3 t"'"(k) = 2t"'"[cos(A;a;) — cos(A;j;)]^ + Iq 
with t-L = O.lt and = ■ We assume an in-plane superconducting order 
parameter A^ with equal amplitude and phase in both layers. The self- 
consistent solution of the mean-field equations then leads to a vanishing 
inter-plane gap A"*"*^ = . 

Magnetic excitations in the superconducting phase are described by 
quasi particles (the fermions) in a BCS-type d-wave pairing state. These 
propagate with effective hopping parameters t , t' strongly reduced from the 
bare parameters t, t' by the small Gutzwiller factor x . The bubble Eq. (Q) 
describes spin-fiip ph-excitations of these particles, which are subject to the 
mode-dependent final-state interaction in Eq. (^), 

J±(q) = aJ(q) ± , J(q) = 2J[cos(5^.) + cos{qy)] (5) 

The inter-plane exchange is chosen as = 0.2 J . The destruction of the 
antiferromagnetic (AF) state of the 1/2-filled system by hole dopingSllUllIll 
is missing in mean-field theory. The necessary correlations of fermions and 
bosons are not contained, and the AF order vanishesat an unphysically 
high doping levelB3 x'^ ~ 0.22 . Therefore we assumeliHlcil a renormalization 
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J aJ of the in-plane nearest-neighbor exchange. Using a = 0.35 re- 
duces down to Xc ~ 0.03 , which is consistent with experiment and makes 
the study of underdoped systems possible. Note that the above-mentioned 
renormalization t ^ t of the quasi-particles comes out of the self-consistent 
calculation, whereas J aJ is a phenomenological model. Our assumn= 
tion of a being independent of doping leads to an AF correlation lengthll3 
^af{x) ~ 1 / A /t — •7^ at T — > , which agrees with known experimentalU 
and theoreticaS results. 



3. RESULTS 

From the susceptibility Eq. (^) the local magnetic excitation spectrum 
is determined from 

ImxfD(w) = JJ Imx=^(q,u;) (6) 

— 7r 

Fig. H shows Imx2D for the odd and even mode at T ^ in the supercon- 
ducting state. A resonance is clearly visible in the odd (— ) mode, at an 
energy uJres = 0.42 J ~ 50meV for x = 0.12 near optimal doping. When x is 
reduced (underdoping) the resonance moves to lower energies and gains spec- 
tral weight. The resonance appears at the same energyllHI as in Imx~(q, 
for fixed wave vector q = {tt,tt) . In addition, both modes lmx2£,{uj) show 
a broad peak ('hump') at energies (jjj^^^p ~ '^hump above ujres ■ In contrast 
to ojres the hump-maxima uj^^^p are almost independent of doping, located 
somewhat below 2A° (2A° « 0.7 J for x = 0.12) . 

The resonance emerges from a pole in Eq. at wave-vector (vr, vr) 
and energy cOres , driven by the effective interaction Eq. (|5|) in the odd (— ) 
mode. Since cUres is slightly below the threshold Oq to the particle-hole (ph) 
continuum, the resonance appears undamped, i.e., as a (5-function. Due to 
the inter-layer coupling J"*" the interaction Eq. (^) is weaker in the even 
(+) mode, and the resonance in Im^"^ is shifted up into the ph-continuum, 
becoming almost suppressed. Consequently, in wave-vector space a sharp 
peak around (vr, vr) is visible only in the odd (— ) mode. This is demonstrated 
in the top panel of Fig. ^. The bottom panel of that figure shows the 
magnetic response at a higher energy close to the hump-maxima oj^^^mp ■ '^^^ 
intensity is much reduced compared to the resonance. However, the magnetic 
excitations at ~ ^hump occupy almost the whole 2D Brillouin zone, and 
despite their small amplitude they contribute to the wave vector integrated 
susceptibility Eq. (^ . The 'hump' can be traced back to ph-excitations 
across the maximum gap A'^ : At q = (vr, vr) the irreducible particle-hole 
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Fig. 1. Wave- vector q integrated odd- and even-mode susceptibilities Imx2D 
for hole filling x = 0.06... 0.18. Parameters are t = 2J , t' = — 0.45t , 
t"*" = O.lt, J"*" = 0.2 J . Curves are calculated with a damping FWHM= 
0.04J ~ 5meV. Inset: q-integrated bubble spectrum Imx^'^ for x = 0.08. 
The maximum is located close to 2A° = 0.78 J . 



bubble Imx*^''"(q, uj) shows a log-type van Hove singularity (vHs) at to = 2A^ , 
remnant of the density of states of the d-wave superconductor. Moving off 
(vr, vr) this vHs splits into three vHs that disperse very weakly throughout 
the Brillouin zone, leading to a soft maximum ('hump') at 2A° in the q- 
integrated (local) bubble spectrum lmx2£){uj) . This is shown in the inset 
of Fig. When the final-state interaction Eq. (|5|) is taken into account, 
the resonance is obtained in the odd mode, and the hump is pulled down to 
'^hump ^ ^hump ^ '^^^ ' ^Iso is the doping dependence of A*^ compensated: 
are independent of doping, while A'' increases with underdoping. 
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Fig. 2. Magnetic response Imx^^^ (in arbitrary units) in wave-vector q 
space for x = 0.08 . Qx , Qy are measured in units of 27r = Ir.l.u. Parameters 
as in Fig. |l]. Top: At the energy uj = u>res where the resonance appears 
in the odd mode. Bottom: At an energy ui close to the 'hump'-maxima in 
both modes. (Note the different amphtude scale.) 



4. COMPARISON TO EXPERIMENT 

Two experimental groups studied the wave-vector integrated magnetic 
response Imxfi) in underdoped YBCO . Refs. 11,2^ reported a line shape for 



YBCOe.e ; which agrees quite well with our theoretical result for x < 0.08 . A 
'hump' in lmx2D (even) appears at « 100 meV , lmx2D (odd) shows a similar 
structure at a somewhat lower energy ~ 90 meV . The well-known resonance 



appears only in lmx2o 34meV. In Refs. ^J2^ two underdoped samples 
YBCOg.? and YBCOg.s have been studied. In the even ("optical") mode 
of YBCOg.? a hump appears around 70meV, whereas the odd ("acoustic") 
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mode shows a weak hump-like structure at ~ 55 meV , separated from the 
resonance at 33 meV . These features tend to move up in energy in the more 
underdoped sample YBCOe.s , while the resonance in lmx2D shifts down to 
25 meV . 

Although the detailed experimental line shapes are not unique, the qual- 
itative features of our calculation are found in the neutron-scattering spectra. 
In particular we reproduce the different dependencies on doping level of the 
resonance at uJres in the odd mode and the hump-like feature at w^^^^ in 
both modes. Also is w^^^p of the odd mode lower than the uj'^ump 
even mode. Theory and experiments can also be compared quantitatively. 
The measured neutron-scattering intensitiesi'0 are of the same order as the 
theoretical ones in Fig. || (using J = 120meV, i.e., Iji^/J = 8.3//^/eV). 
The maximum of the hump in the even, odd mode in Fig. |l| occurs at 
0;+'^ 0.6 J, 0.53 J = 72meV, 64meV, in good agreement with the mea- 
surements Ref. |9| on YBCOq.t at low temperature. 




hole filling x 

Fig. 3. 2|;U-^| , 2A^ , the ph-threshold VLq , and the resonance energy u^es as 
function of doping at T ^ . 



5. CONCLUSION 

The local magnetic excitation spectrum Inrx^oi'^) characterized by 
two energy scales that behave differently with hole filling x (doping). The 
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resonance energy tUres follows the particle-hole threshold ^Iq , which in un- 
derdoped systems is determined by the chemical potential //-^ of the fermions 
as OJxes < = 2|//-^| . The maxima of the humps, on the other hand, are 
determined by the gap A*^ through w^^^p ~ 2A'^ . When x is reduced from 
optimal doping, and therefore ujres decrease quickly, while A*^ increases. 
This is displayed in Fig. ^ . It has been noted above that the mean-field 
theory describes magnetic excitations in terms of quasi particles (QP) (the 
fermions) with a reduced Fermi velocity vp ^ {x + 0.15J/t)vF . Hence in 
underdoped systems the QP's chemical potential comes out (much) smaller 
than the gap, < A*^ , and thus determines the scale for u>res ■ This leads 
to the observed decoupling of the resonance energy from the gap A^ , while 
A*^ is still visible through the 'humps' in the local spectrum lmx2£) ■ 
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